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APPLICATION  OF  STATISTICAL  ESTIMATION  TECHNIQUES  TO 
GROUND-BASED  PASSIVE  PROBING  OF  THE  TROPOSPHERIC 
TEMPERATURE  STRUCTURE 


by 

Ed  R.  Westwater 
and 

Otto  Neall  Strand 
Abstract 

The  theory  of  radiative  transfer  in  the  microwave  region  is 
discussed.  A  method  foF*fche  numerical  solution  of  a  Fredholm  integral 
equation  of  the  first  kind  is  derived  and  illustrated.  The  method 
employs  an  a  priori  constraint  vector  together  with  covariances  of 
both  the  constraint  vector  and  the  measurement  enors.  The  method 
automatically  incorporates  an  optimum  amount  of  smoothing  in  the 
sense  of  maximum -likelihood  estimation.  The  problem  of  obtaining 
optimum  basis  vectors  is  discussed.  The  trace  of  the  covariance 
matrix  of  the  error  in  the  solution  is  used  to  estimate  the  accuracy  of 
the  results.  This  trace  is  used  to  derive  a  quality  criterion  for  a 
set  of  measurements  and  a  given  set  of  constraint  statistics.  Examples 
are  given  in  which  the  behavior  of  the  solution  as  obtained  from  a 
specific  integral  equation  is  studied  by  the  use  of  random  input  errors 
to  simulate  measurement  errors  and  statistical  sampling.  The  quality 
criterion  and  behavior  of  the  trace  of  the  error  covariance  matrix  for 
various  bases  are  also  illustrated  for  the  equation  being  considered. 

A  least-squares  analysis  of  microwave  absorption  data  to  determine 
line  width  for  water  vapor  and  oxygen  is  presented.  Calculations  of 
the  error  covariance  matrix  and  quality-criterion  values  to  be  expected 
when  probing  for  the  tropospheric  temperature  structure  from  micro- 
wave  emission  measurements  are  presented. 


1.  Introduction 


The  general  subject  of  passive  probing  of  the  atmosphere  by 
electromagnetic  means  has  recently  received  attention  by  several 
authors:  Fleming  and  Wark,  [l]  ;  Kaplan  [  2]  ;  and  King,  [3]  . 

In  spite  of  the  wide  differences  in  the  specific  techniques  applied 
to  each  spectral  region,  several  elements  are  common  to  all  pas¬ 
sive  probing  schemes.  The  first  of  these  elements  is  the  sensor, 
which  measures  relevant  properties  of  the  radiation  emanating 
from  the  spatial  region  being  probed.  Secondly,  detailed  know¬ 
ledge  is  required  concerning  the  mechanisms  by  which  radiation 
is  generated  at  the  region  being  probed  and  by  which  the  radiation 
is  propagated  to  the  sensor.  The  last,  but  certainly  not  the  least, 
is  the  manner  by  which  information  concerning  the  physical 
property  of  interest  is  derived  from  the  set  of  radiation  measure - 
ments,  the  so-called  "inversion  process". 
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The  probing  methods  discussed  here  attempt  to  derive  informa¬ 
tion  about  the  temperature  distribution  of  the  troposphere  by  measur¬ 
ing  the  microwave  brightness  temperature  of  the  atmosphere  from 
the  ground  looking  up.  The  propagation  and  emission  of  this  radia¬ 
tion  are  determined  for  the  most  part  by  the  absorption  coefficients 
of  oxygen  and  water  vapor. 

The  inversion  method  we  derive  here  is  illustrated  for  the 
physical  situation  of  ground-based  upward  probing;  however,  the  in¬ 
version  equations  are  applicable  to  more  general  situations.  This 
technique  extends  methods  given  by  Phillips  [4],  Twomey  [  5  ] ,  [  6]  , 
and  Twomey  and  Howell  [7]  for  the  solution  of  integral  equations 
of  the  first  kind.  These  methods  use  a  controlled  amount  of  smooth¬ 
ing  in  solving  the  matrix  system  derived  from  a  quadrature  approxi¬ 
mation,  but  give  no  systematic  method  of  determining  the  required 
amount  of  smoothing.  In  this  paper,  we  derive  and  illustrate  a 
general  least-squares  process  for  estimating  the  solution  and  show 
that  it  contains  Twomsy's  method  of  weighting  an  a  priori  constraint 
vector  as  a  special  case.  The  optimum  amount  of  smoothing 
(in  the  sense  of  maximum -likelihood  estimation)  depends  on  certain 
covariance  matrices  describing  the  measurement  process.  In  a 
manner  similar  to  that  of  the  more  recent  work[8J  ,  [9],  [  1 0 J ,  and 
[  1 1 3 ,  we  derive  and  use  statistically  orthogonal  basis  vectors. 
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Certain  other  bases  are  also  used  and  numerical  comparisons  are 
made  for  a  particular  kernel.  In  our  method,  the  covariance  matrix 
of  the  error  in  the  resulting  solution  is  calculated  and  its  trace  is 
used  as  a  measure  of  the  error  in  the  final  results.  It  is  shown  in 
this  paper  that  the  trace  of  the  error  covariance  matrix  is  mini¬ 
mized  when  the  number  of  basis  vectors  equals  the  number  of 
quadrature  points.  We  use  the  resulting  minimum  trace  as  a  quality 
criterion  describing  the  effectiveness  of  a  given  matrix  equation 
(with  specified  accuracy  of  observation  and  a  specified  distribution  of 
the  constraint  vector)  in  reducing  the  variance  of  the  a_  priori  con¬ 
straint  vector.  We  also  derive  and  discuss  the  condition  for  an 
optimum  basis  and  present  numerical  examples  using  random  input 
errors  to  simulate  the  measurement  process.  We  discuss  and 
present  a  least -squares  analysis  of  line  width  and  absorption  data 
for  both  water  vapor  and  oxygen.  We  also  present  numerous  calcu¬ 
lations  of  the  solution  quality  criterion  and  indicate  how  this  criterion 
could  be  used  to  determine  both  experimental  feasibility  and  optimum 
choice  of  probing  frequencies. 


Z.  Theory  of  Radiative  Transfer  in  the  Microwave  Region 

The  equation  of  transfer  in  a  steady  field  of  radiation  in  the  case 
of  local  thermodynamic  equilibrium  is  given  by  Kondrat’yev  [12] 
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where  1^  is  the  intensity  of  radiation,  Bis  the  zenith  angle,  h  is  the 

height  above  the  surface,*  and  c  are  the  mass  absorption  and  mass 

V  V 

scattering  coefficients  respectively,  p  is  the  density  of  material  absorb¬ 
ing  the  radiation,  B^is  the  Planck  intensity  function,  and 

p^  (h,  r'  ,  r)  is  the  fraction  of  energy  incident  upon  the  scattering 
volume  in  the  direction  r1  which  is  scattered  in  the  direction  r. 

The  integration  over  dui'  extends  over  all  directions  of  r'.  For 
a  cloudless  atmosphere  and  for  microwave  frequencies,  the  scatter¬ 
ing  coefficient  is  zero  and  the  equation  of  transfer  for  a  purely 
absorbing  atmosphere  is  obtained 
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for  the  intensity  received  "looking  up"  at  a  height  h  and  at  a  zenith 
angle  9.  In  (3),  1^ (°° )  is  the  unattenuated  intensity  incident  outside 
the  atmosphere  and  t  (h,  ®)is  the  total  vertical  attenuation  given  by 


5 


(4) 


CO 

i,«)  =  ^KvPdh' 


The  Planck  function,  B  ,  multiplied  by  the  absorption  coefficient  is  pro¬ 
portional  to  the  power  flux  density  emitted  per  unit  bandwidth  and  may  be 
approximated  by  the  Rayleigh -Jeans  law  in  the  microwave  region  [  13]  . 

B  ^kl,  (5) 

^  v  2 

c 

where  k  is  the  Boltzman  constant,*:  is  the  speed  of  light,  and  T  is  the 
absolute  temperature.  The  intensity,  1^  may  be  related  to  the  bright¬ 
ness  temperature,  T^  ,  by  the  defining  relation 


I  -^kT. 

V  c2  v 


It  follows  immediately  from  (3),  (5),  and  (6)  that 


Tb  (h)  =  TbC)  e 


-  TJh,a)/cos  9  ^ 


-  W  p  dh"/cos  9 

J  v 


+  jKvpTe 


dh‘ 
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The  apparent  sky  temperature,  as  measured  by  an  antenna 
immersed  in  a  radiation  field,  will  be  the  contribution  of  (7)  from 
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each  direction  seen  by  the  antenna  multiplied  by  the  appropriate 
antenna  gain  factors  and  integrated  over  the  antenna  pattern. 

If  it  is  desired  to  determine  the  brightness  temperature  of  the 
atmosphere  from  directions  confined  to  a  narrow  solid  angle,  it  is 
necessary  to  have  a  highly  directional  antenna  with  low  side  lobes 
such  that  the  contributions  from  the  desired  direction  are  maximized. 
The  noise  temperature  of  the  antenna  can  be  measured  directly  by 
comparing  the  received  sights  1  with  that  of  a  known  source. 

In  using  (7)  to  deduce  T(h)  from  measurements  of  T^»  one  may 
vary  either  zenith  angle,  G,or  frequency,  v.  In  ground-based  probing, 
however,  there  are  reasons  to  avoid  the'angular  scheme.  First, 
an  assumption  of  horizontal  homogeneity  of  the  atmosphere  is  necessary. 
It  is  also  probable  that  some  of  the  space  swept  out  by  the  observing 
antenna  will  include  discrete  cloud  structures  whose  emission  is  not 
easily  taken  into  account.  Finally,  it  is  very  difficult  to  take  account 
of  ground  emissions  and  reflections  because  of  uncertainties  in  antenna 
gain  patterns.  We  will  therefore  only  consider  inversion  using  measure¬ 
ments  at  different  frequencies.  This  scheme  reduces,  but  does  not 
eliminate,  the  above  difficulties.  Henceforth,  it  will  be  assumed  that 
no  clouds  are  present  in  the  main  beam  of  the  observing  antenna. 
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Equation  (7)  is  non-linear  in  the  temperature  because  the 
absorption  coefficient  depends  o.i  temperature  and  pressure. 

As  discussed  by  Westwater  [  14],  this  non-linearity  can  be  handled  by 
assuming  a  mean  temperature  profile  to  determine  a  first  approxi¬ 
mation  to  the  absorption,  solving  the  resulting  Fredholm  equation  of 
the  first  kind,  and  iterating.  The  next  section  presents  a  method  of 
solving  the  linearized  equation  in  the  presence  of  measurement  error 
by  statistically  weightings  priori  knowledge  of  the  temperature  pro¬ 
file  against  information  given  by  the  integral  equation.  Because  of  the 
nature  of  the  constraint,  it  is  presumed  that  convergence  of  the  itera¬ 
tion  sequence  will  cause  no  problem,  $ven  in  the  presence  of  large 
experimental  errors.  This  has  not  yet  been  verified. 

3.  Inversion  Method 

Let  {f(y),  a  ^  y  ^  b}  be  a  Gaussian  random  process  having 
continuous  sample  functions.  If  K(x,  y)  is  a  continuous  function  of 
x  and  y  and  if 

b 

J  K{x,  y)  f  (y)  dy  =  g  (x)  (8) 

a 
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then  g{x)  is  also  a  Gaussian  process  with  continuous  sample  functions. 

The  main  idea  is  to  estimate  f(y)  by  measuring  g(x)  at  various  values 

of  x,  say  x  =  x.,  i=  1,  .  .  .  n.  Introducing  a  quadrature  formula  of  the 
I P  xm 

type  J  h(y)dy=  £  w.h(y.)  gives  a  matrix  equation 
a  J--1  J  J 

Af  =  g  (9) 

where 


A  =  ( A„ )  i  =  1 , 2 ,  ...  n;  j  =  1 , 2 ,  ...m, 
m  is  the  number  of  quadrature  abscissas, 
n  is  the  number  of  observations  of  g(x), 
y^  are  the  quadrature  abscissas, 

x.  are  the  specific  values  of  x  for  which  g(x)  is  observed, 
w.  are  the  weights  associated  with  y., 


f. 

J 


gi  =  gfcj). 

T 

f  =  [f,f_,  .  .  .  f  ]  is  the  column  vector  of  unknown  functional  values, 

12  mJ 

(the  superscript  T  denotes  matrix  transposition  throughout  this  paper), 

T 

and  g  =  [g,g0...g]  is  the  column  vector  of  values  of  g(x). 

1  2  n 

We  assume  that  the  mean  vector  E(f)  =  f  and  the  covariance  matrix, 

o 

Sjj,  of  f  are  known.  (E(  )  denotes  the  expected -value  operator  and  S 
denotes  the  covariance  matrix  of  any  vector  v  throughout  this  paper). 


By  the  linearity  of  E  and  the  propagation  rule  for  covariance  matrices. 

T 

respectively,  we  have  E(g)  =  g  =  Af  and  S  =  AS.A  .  In  the  cases 

o  o  g  f 
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of  interest  here  g  is  a  vector  of  measurements  subject  to  error. 


Thus  one  observes 


ge  =  g  +  e 


instead  of  g.  In  addition  to  the  assumptions  already  made,  we  assume: 

a)  the  errors  are  independent  of  f,  hence  independent  of  g; 

b)  the  errors  e.  are  normally  distributed  with  zero  mean  and 

known  covariance  matrix  S  ; 

e 

c)  the  quadrature  errors  are  negligible  with  respect  to  e  ; 

d)  S  and  S  are  both  nonsingular  with  dimensions  nxn  and 

£  f 

mxm,  respectively. 

It  is  convenient  in  practice  to  represent  the  solution  in  terms 
of  a  small  number  of  basis  vectors.  We  let  U  =  [  |  |...  1 1^] 

(lines  denote  matrix  partitions  throughout  this  paper)  be  a  matrix 
of  k  linearly  independent  m -component  basis  column  vectors  and 


c  =  [Clc2...  C^J 


be  a  vector  of  coefficients  to  be  determined 


where  k  ^  m. 

The  inference  problem  to  be  solved  is  the  following.  Given  the 

observed  vector  g  and  the  covariance  matrices  S,  and  S„  ,  estimate 
°e  f  e 

f  and  c  under  the  constraint  gg  -  Af  =  e.  Putting 

ri  =  f-fo,  h  =  g  -  gQ  -  A(f  *  AT! 


h  =  h  +  s  =  AT]  +  s, 
e 
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where  r|  is  to  be  approximated  by  UC  gives  an  equivalent  problem: 

given  hg  ,  estimate  T]  and  e  in  the  form  subject  to  the  constraint 

h  -  AUC  =  €  .  First,  note  that  the  maximum -likelihood  estimate  of  T| 
e 

in  the  absence  of  measurements  g^  is  r|  =  0.  The  estimate  T)  =  0  may 
be  coupled  with  the  set  of  measurements,  g^,  to  form  an  "effective" 

r  o  i 

measurement  vector  jj- — j.  From  this  vector  and  the  Gauss -Markov 
least -squares  principle  for  correlated  variables  [15]  ,  a  maximum- 
likelihood  estimate  of  I  —  j  can  be  obtained.  Since  the  covariance 

matrix  of  is  |~  ^  -|  ^ —  ,  the  desired  estimate  j  ~|p*J  must  be  such 

e 

that  R(C,  e)  is  rendered  stationary,  where 


F.(C,  e)  s  (UC)T  S  '1  (UC)  +  eT  S e  +  pT(AUC  -  h  +  e),  (12) 

X  G  C 


and  (j.  is  a  row  matrix  of  Lagrange  multipliers.  Differentiating  with 
respect  to  C,  Cjand  p  and  equating  to  zero  gives  the  solution 


where 


UC  =  UD_1  UTAa  s  _1h 

e  e 


S  =  h  -AUD_1UTATS  _1h  ? 

e  e  e 


D  =  UTXU  and  X  =  S  _1  +  ATS  _1A.  (15) 

X  c 
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We  may  also  rewrite  (13)  in  the  equivalent  form 

UC  =  U{UTSf'1U)  _1  UTAT[s,*  AU(UTSf"2  U)'1UTAT]'1he 

We  apply  the  propagation  rule  for  covariance  matrices  to  UC-7 
where  UC  is  given  by  (13): 

uc  -  r*  =  ud'1  utat  S  -1(An  -  £  )  -  ti 

€ 


and 


=  [UD_1  UTATS£_1  A-l]sf  [UD_1  UTATS„  ^-Ij1 


UC  -  n 


-1_ 


-1  ,  T  0  *  11  T 


[UD_1  UTATS£  jSs  j_UD  *  U*A*  S£ 


f. 


so  that 
S 


UC 


_  TYS  Yt  -  (TYSf  S£  YT)  x  Sf  TYT 


where 


T  -1  -IT 

Y  =  A  S  A  and  T  =  UD  U  . 
e 


-1 


Substituting  Y  =  X  -  Sf  into  (17)  and  noting  that  TXi  -  T  give; 


s  Sf  +  T+T  (XS£X)  T  -  TXSf  -  SfXT, 


and 


S-  =  AS. .  _  A 

€  -  €  UC  -  n 


.  (16) 


(17) 


(18) 


(19) 


(20) 
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Expressions  (13),  or  (16),  and  (14)  give  the  desired  statistical  estimates 

of  the  solution  for  any  given  basis  U.  Equations  (19)  and  (20)  give  the 

error  covariance  matrices  by  which  the  quality  of  the  results  may  be 

estimated,  as  will  be  seen  in  the  next  section.  In  the  work  reported 

in  this  paper  no  use  has  been  made  of  the  estimate  (14)  of  e  of  the 

expression  (20)  giving  S* 

€  -  0 

Various  special  ca^jgs  may  be  noted.  If  k  =  m,  that  is,  if  U  is 
a  nonsingular  mxm  matrix,  (13),  (16),  (19),  and  (20)  become, 
respectively, 


A  -1  T  -1 

UC  =r\=x  1  As  h  , 
e  e 

(21) 

*1  =  S  AT  H  " 1  h 
f  e 

(22) 

=  x"1 

(23) 

-1  T 

Sa  S  AX  a 

e  -  e  ’ 

(24) 

where 

H  =  S  +  AS  AT  . 
s  f 

From  (22),  a 

result  equivalent  to  (23)  is 

x"1  =  Sa- 

r,-r\ 

=  SfPSfPSf  -  2SfPSf  +  Sf  +  Sf  ATH  _1Se  H_1ASf  , 

(25) 

where 

T  -1 

P  =  A  H  A. 
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The  study  of  special  cases  will  be  confined  to  the  case  k=m  here.  If 
S  “•  0,  i.  e.  ,  if  the  measurements  have  large  errors,  we  obtain 

A  _1 

11  =  0  or  f  =  f  and  Sa  =  S,.  If  S,  -*  0,  then  X  will  be  nonsingular 
o  T)  -  T|  f  f  ° 

A  T  -1  -IT  -1 

if  and  only  if  rank  A=m.  Then  r|  =  (A  S  A)  AS  h  and 

see 

T  -1  -1 

S*  =  (A  S„  A)  .  If  S  -*0,  i.e.,  if  the  measurements  of  g  are 
T|  - 1|  e  e  6 

perfect,  and  if  rank  AU  =  n  then  n^m  and  (22)  and  (23)  become 

respectively,  r\=  SfAT"^ASf AT  j  1  h  and  S*  =  Sf  -  SfAT(ASfAT)  _1ASf. 

In  computational  practice  the  latter  situation  will  not  occur  since  the 

quadrature  errors  will  invalidate  the  assumptions  as  S  “O.  Finally, 

if  S^-*0,  then  (22)  and  (25)  give,  respectively,  *  =  0  or  f  =  f  and 

S-  =  0. 

’I  -  ^ 

We  may  show  the  correspondence  of  the  present  method  with 

7  2 

that  of  Twomey  [6]  as  follows:  for  k=  m,  let  S  =  a  “  I,  S  =a  I, 

€  €  xi 

2/2  _  a  — 

V  =  ae/  °£  '  =  P«  +T1=  f;and  gg  =  g  in  (21)  to  correspond  with 

Twomey' s  notation  [  6]  .  Then 


f  T 

f  =  j_A  A  + 


v-rt 


T_ 

A  g  +  v  p 


(26) 


which  is  identical  to  Twomey's  equation  (6),  p.  105.  Thus,  if  the 
covariance  matricas  Sf  and  Sg  are  both  scalar,  the  optimum  choice 
of  y  (in  the  sense  of  maximum -likelihood  estimation)  is  given  by  the 
ratio  of  variances.  However,  in  commonly  occuring  physical  situations 
has  substantial  off-diagonal  elements. 


The  basis  U  =  I  with  m  =  k  gives  an  adequate  solution  via  (21), 

(23),  and  (24).  However,  it  is  often  convenient  to  represent  the 
solution  in  terms  of  k  parameters,  where  k  <  m,  using  (13),  (19),  and 
(20).  The  quality  of  such  a  representation  depends  on  the  choice  of 
basis.  In  the  next  section  we  discuss  the  problem  of  basis  optimi¬ 
zation. 


4.  Determination  of  Optimum  Basis  Vectors 

Previous  solutions  by  Wark  and  Fleming  [  10]  ,  and  Twomey  [  11] 
have  employed  basis  vectors  (i.  e.3columns  of  U)  which  were  obtained 
to  fit  the  data  in  an  optimum  manner  without  reference  to  the  integral 
equation  to  be  solved.  It  will  be  shown  below  that  the  integral  equation 
(as  represented  by  the  matrix  A  which  is  an  ingredient  of  ) 

plays  a  part  in  obtaining  the  optimum  matrix  U,  but  (see  section  5) 
the  use  of  an  exactly  optimum  basis  U  is  often  not  critical,  so  that 
U  as  employed  in  earlier  solutions  should  suffice.  This  classical 
determination  of  U  (as  exemplified  by  Obukhov's  paper  [16],  for 
example)  is  presented  here  and  interpreted  in  terms  of  the  covariance 
matrix  of  the  residuals.  The  criterion  which  must  be  satisfied  for  an 
optimum  basis  for  the  integral  equation  is  then  easily  seen. 


Consider  the  quadratic  form 


(d  -  UC)  (d  -  UC)=  q  (U,  C). 


(27) 
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Here  U  and  C  have  the  same  meaning  as  in  section  3  and  d  is  a 


sample  vector  from  a  population  having  a  multivariate  m  -dimensional 

normal  distribution,  say,  with  zero  mean  and  covaariance  matrix  S,. 

d 

For  any  specified  basis  U,  q(U,  C)  is  minimized  when 

UC  =  Ull^U)'1  Ur  d,  (28) 

as  can  be  seen  by  differentiation.  Thus  the  minimum  value  of  q(U,  C) 
for  any  gi\en  basis  U  is  given  by 

T  f  T  -IT  “IT”  T  -IT  1 

q  .  (in  =  d  U(U  U)  U  -I  I  |U(U  U)  U  -lid.  (29) 
min.  L  J  U  J 

Various  criteria  may  be  used  to  determine  what  is  meant  by 

"smallness"  of  q  In  the  classical  case  the  optimum  U  is  that  which 

minimizes  the  expected  value  of  q  (U).  Let  TrB  =  trace  B  denote 

main 

the  sum  of  the  diagonal  elements  of  any  matrix  B.  By  expanding 
and  employing  the  linearity  of  the  operator  E  we  may  show  that 

Ef  dXQd]  =  Tr  Q  (30) 

a 

for  any  real  symmetric  matrix  Q.  Applying  (30)  to  (29)  gives  U 
as  that  basis  which  minimizes 

Tr{Sdl  I  -  U(UTU)_1U^j  } 

=  T rS  -  Tr  UTSJU(UTU)"i  .  (31) 

a  a 
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In  the  last  expression  we  have  used  the  relation  Tr  AB  =  Tr  BA,  which 

holds  whenever  AB  is  square.  Expression  (31)  is  invariant  with  respect 

to  replacement  of  U  by  UB,  where  B  is  any  nonsingular  kxk  matrix. 

T  -IT 

In  fact,  the  matrix  S,U(U  U)  U  is  unaffected  by  such  a  transforma- 

d 

tion.  Furthermore,  let  J  be  an  mxm  symmetric  positive  definite 

T 

matrix  and  let  V  be  any  mxk  basis  matrix.  Since  V  JV  is  positive 
definite,  there  exists  [  If*]  a  real  nonsingular  upper  triangular  kxk 
matrix  B  ^  such  that  V^JV  =  (3  ^)^(B  *)•  Let  U  =  VB.  Then 
UTJU  =  BT(VTJV)  B  =  I.  Thus  any  normalization  of  the  form  LTAU  =  I 
may  be  assumed  without  loss  of  generality.  We  suppose 


U  U  =  I. 


(32) 


T 

Hence  U  must  be  found  which  maximizes  Tr  (U  S,U}  subject  to  the 

d 

constraint  (32).  By  partitioning  U  into  individual  columns  of  basis 
vectors  and  considering  the  ma.ximization  process  involved,  it  can 
be  shown  [  18]  that  the  columns  of  U  must  be  chosen  as  those 
eigenvectors  corresponding  to  maximal  eigenvalues  X  a  .  .  .  a  X.  >  0 

X  iC 

of  S^,.  Thus  U  must  satisfy 


where 


SJ  U  =  U  A 
d 


(33) 


A  =  diag  [ 
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is  sometimes  called  [8j  the  fraction  of  the  total  variance  "explained" 
by  the  basis  U. 

It  follows  from  (28)  and  (31)  by  us.  the  propagation  rule  for 
covariances,  with  TrAB=  TrBA,  that 


E  £q  .  (U)  }  =  T rS  .. 

nun  UC  -  d 

Thus  E  {q  .  (U)  3  may  also  be  given  another  interpretation.  As  a 
^nun  '  r 

criterion  of  the  size  of  the  errors  resulting  from  a  distribution 

having  zero  mean  and  second  moments  defined  by  a  given  covariance 

matrix  S,  we  may  choose  TrS.  Since  the  trace  is  equal  to  the  sum 

of  the  eigenvalues  of  a  matrix,  this  interpretation  is  geometrically 

equivalent  to  defining  the  "size"  of  an  m  -dim ensional  ellipsoid  as  the 

sum  of  the  squares  of  its  semi-axes.  The  interpretation  of 

E{q  •  (U)  ]  as  the  trace  of  the  covariance  matrix  will  be  freely  used 
nrun 

throughout  the  remainder  of  the  report.  Note  that  E^q  (U)  3  is  also 

'^min 

equal  to  m  times  the  expected  overall  mean-  square  error  of  the 
approximation  d  —  UC. 
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Although  no  use  is  made  of  the  result  in  what  follows,  it  is  interest¬ 
ing  that  the  minimum  of  the  quadratic  form 

(U,C)  -  (d  -  UC)rs/  L(d  -  UC) 

d  d 

has  the  expected  value  m  -k  for  any  choice  of  U,  and  therefore  cannot 
be  used  to  determine  an  optimum  basis. 

From  the  preceding  discussion  it  follows  that  an  optimum  basis 
U  for  use  with  the  integral  equation  would  be  that  which  minimizes 

TrSuc_n  *  Tr{Sf  +  T+T  (XS£X)  T-TXSf-SfXT  J  ,  (35) 

From  the  definition  ot  T  we  see  that,  even  with  the  normalization 
T 

U  XU  =  I,  expression  (35)  is  of  fourth  degree  in  the  elements  of  U. 

Since  minimization  of  (35)  is  difficult,  and  since  classical  basis 
vectors  seem  to  sulfice  in  practice,  the  determination  of  U  from 
(35)  will  not  be  pursued  further  in  this  paper. 

5.  Quality  Criterion 

The  quantity  TrX  ^  can  be  used  as  a  measure  of  the  accuracy  to  be 

expected  from  a  given  integral  equation  with  given  covariance  matrices 

S,  and  S  and  a  given  set  {x.  }  of  values  of  x.  For  k  =  m,  i.  e.,for  a  non- 

singular  mxm  matrix  U,  S  reduces  to  X  ^  by  (23).  Furthermore,  as 

demonstrated  below,  TrS,._  STrX  ^  for  any  basis  U.  Thus,  TrX 

U\J  “  M 

is  a  measure  of  the  best  that  can  be  done  for  a  given  problem.  Calcu¬ 
lating  TrX  *  and  comparing  with  TrS£  indicates  the  amount  of  improve¬ 
ment  over  the  a  priori  statistical  knowledge  that  can  be  expected  in  a 
given  case. 
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£7 

K 


*VrgpB>  4* 


-l 


Th  eore’ii:  Tr  S.|r  .  -TrX  s  0.  where,  as  previously  defined, 

r  ■  U  Vj  ~  ’ 


-1  T  - 1 

x  =  s.  +als  a, 
1  £ 


D  =  U  XU, 

-1  T 

T  =  UD  U  , 


Tr  S 


uc  _  -  TrX"1  =  Tr  [Sf  T  +  TXSfXT  -  T  X  S£  -  SfXT  -  X_1  j 


U  = 


[ui|ua|---  |  uk] 


U  ...  U  m  where  the  U.  are  linearly  independent 
2  k  l 


and  m  s  k. 

T  T  -1 

Lemma:  Tr  B  S£  B  .  TrB  X  B  where  B  is  any  real  mx  m  matrix. 
Proof:  Let  S£=(Sj.),  X  *  -  (X.  ^),  and  B  =  (B..).  By  direct  calculation, 


m  m  m  «  m 

..B  =  7  B  S  B  .  »  V  (J  T  Sp,  B  B  ]  =VbTS.  B. 

1  i,Piq  pi  p«.  i‘  iVi\P4i  q4ip<1  qs  py  t4i‘  ■  1 


T 

Tr  B  S 


w 


here  B  is  the  i— —  column  vector  of  B.  By  an  extension  [  19 ]  of 


Bergstrom's  inequality  we  have,  for  any  real  column  vectors  y  and  z, 

T  T  - 1  ,  T  2 

(y  Sf  y)(z  Sf  z  )  -  (y  z  j  . 

-1  T-l  T  T  - 1 

Since  X=  S.  +  A  S  A  so  that  2  Xz  ^  z  S  z,  we  have 

♦  r  I 
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(yTSfy}(*TXz)  ^  (yTz)^  Pitting  y  =  B.,  z  =  X  1  B.  gives 

(B  ^  S,.  B  )(B^  X  *  B.)  -  (B  *  X  *  B.)2.  If  B.  P  0,  dividing  out  the  po^iti\ e 
a  f  \  1  i  i  i  i 

quantity  bT  X.  ^  B.  gives  Ik  S^B.  -  B^  X  B^.  (If  =  0,  the  in¬ 
equality  is  trivially  true.)  Therefore, 
m  m 

T  V  T  r1  T  T  -1 

TtB  S.B*  Bi  S,B.  s  Bj  X  B.  =  TrB  X  B.  This  proves 
f  .  ,  f  1  T  1  1 

1=  1  1=1 

the  lemma.  It  can  now  be  applied  to  prove  the  theorem.  Note  that 
Sue  q  "X  1  may  be  written  as  ^  -  X  1  =  (XT  -  I)  TSf(XT  -  I  )  + 

T  -X  By  the  lemma 

Tr  [S  T  _  -  X*1]  a  Tr  {(XT -I)TX_1  (XT  - 1)  j-TrT-TrX'1 
=  Tr  [TXT  -  2T  +  X-1  -  T  -X"1  }.  Since  TXT  =  T  we  have 

Trisuc-n  T'-d- 

There  are  conditions  for  which  a  basis  may  be  chosen  such  that 

Tr  ST_£  ^decreases  as  k  increases,  and  approaches  TrX  ^  for  k 

considerably  less  than  m.  Such  behavior  is  illustrated  in  f  igure  1, 

which  presents  the  values  ox  TrSn  vs.  k  for  the  integral  equation 

(36),  given  and  discussed  in  section  6-  The  conditions  are: 

S  =  .  01  L  n  -  7,  m  =  15  with  {x.}  and  £y.}  as  given  in  section 
e  i  J 

6*  The  Obukhov  basis  resulting  irom  the  Sf  0f  table  I  is 
used.  Figure  2  oreser.ts  the  values  of 
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TrSjjc  vs.  k  under  similar  conditions  using  the  power  basis 

U..  r  y"!  and  figure  3  presents  similar  results  for  the  trigonometric 

basis  U.  -  sin  ^  i  =  1,  2,  ...  15,  j  =  1, 2,  .  7 

Ui,  2j  +  1  *  COB  CbT^a  )  1  =  l’  Z>  *  ’  ‘  15,  j  =  °*  l’  2*  •’  ‘  '  7' 

Evidently  k  =  6  will  suffice  for  the  power  basis  and  k  =  11  will  suffice 

for  the  trigonometric  basis.  Figures  1,2,  and  3  are  combined  for 

easy  comparison  in  figure  4. 

The  results  of  figures  1  through  4  arise  from  reasonably  well- 
chosen  bases.  The  type  of  behavior  to  be  expected  for  ill-chosen  bases 
is  indicated  in  figures  5  and  6.  Figure  5  presents  TrS^^  ^vs.  kiora 
basis  consisting  of  random  numbers  uniformly  distributed  between  -1 
and  +1.  Figure  6  presents  similar  results  for  the  basis  U„  =  Clearly, 
in  these  cases  nothing  less  than  a  full  complement  of  basis  vectors 
(i.  e.,k  =  m  =  15)  will  achieve  the  desired  accuracy.  By  noting  that 

TrS  is  m  times  the  theoretical  overall  mean-square  error  of  the 

UO  ”  -I 

fit,  we  see  that  TrS^.  -  120.407  corresponds  to  an  overall  rms  error  of 
2.  88  and  that  TrX  *  =  23.  946  corresponds  to  an  overall  rms  error  of 


1.26.  This  shows  the  extent  to  which  the  original  statistics  of  f  may 


be  improved  by  using  the  integral  equation  in  this  particular  case. 
The  analysis  presented  here  serves  a  similar  purpose  to  the 
analysis  of  the  "degree  of  independence"  of  the  measurements  of 
g(x)  as  discussed  by  Twomey  [  6]  ,  [ll]  •  However,  in  the  present 
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case  the  actual  rms  error  to  be  expected  is  found  and  the  results 

are  necessarily  dependent  on  the  statistics  S  and  S..  as  well  as 
*  ct 

on  the  kernel  of  the  integral  equation. 


6.  Numerical  Inversion  Results 

For  purposes  of  numerical  experimentatioi  the  following  equation 


was  used: 


-aWHo  aw  je  'y/Ho  «„)  „a(x)H„  e’y/H° 


f(y)  e  o 


dy  =  g(x). 


This  equation  occurs  in  remote  atmospheric  probing  work  if  an 
exponential  atmosphere  is  assumed.  Here  H  is  a  constant  and 


o(x)  is  given  by 


a(x)  =  1 .  1  x  -  1 


for  the  purposes  of  this  section.  If  f  =  f  s  const,  is  inserted  in  the 
left  side  of  (36)  and  the  integrations  performed,  the  resulting  right 


side  is: 


g(x)  = 


e-a(x)Ho(l  -e 


-H/Ho>] 


Substituting  (38)  into  (36)  will,  of  course,  give  an  integral  equation 

whose  correct  solution  is  f  =  f^.  In  the  work  reported  here,  the 

values  H  =  10,  H  =  5  and  f  =  293.997  were  assumed.  To  study  the 
o  *  c  7 


~  *gS- ■  wfot frgMHffifrgy 


I.  COVARIANCE  MATRIX  OF  t 
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behavior  o£  (36)  in  the  presence  of  measurement  error,  random 

errors  distributed  in  accordance  with  S  were  added  to  g(x)  as  given 

€ 

by  (38)  and  the  mean  value  f  was  varied  randomly  from  f  in  a  manner 
1  o  c 

determined  by  S^.,  which  is  given  in  Table  1. 

Since  the  solution  f  simulates  an  atmospheric  temperature  pro¬ 
file,  was  obtained  from  radiosonde  data  from  240  soundings  during 
August  at  Denver,  Colorado.  The  matrix  was  assumed  to  be 
scalar  and  the  values  assumed  will  be  indicated.  The  method  of 
constructing  random  errors  is  indicated  below. 

The  CDC  3600  computer  has  a  F  OR  TRAN  function,  RANF  (-1), 
which  gives  random  numbers  uniformly  distributed  between  0  and  1. 
These  numbers  are  random  in  the  same  sense  that  a  table  of  random 
numbers  is  random;  that  is,  if  the  series  of  statements  Y.  =  RANF(-l), 
i=l,2,...M  is  repeated  in  the  same  manner,  where  i=  1  refers  to 
the  first  use  of  RANF  (-1)  in  a  given  computer  run,  the  same  sequence 
will  result,  but  the  sequence  itself  is  random.  It  follows  from  the  cen- 

tral  limit  theorem  [20]  that  the  variable  V,,,  where 

M 


VM/12 


-  M/2 


(39) 


is  approximately  a  random  normal  deviate  (i.e.^it  has  a  normal 
distribution  with  mean  zero  and  standard  deviation  1),  the  degree  of 
the  approximation  depending  on  the  size  of  M.  A  sample  of 

A 

9,  000  values  of  Vj^  was  calculated  for  M=l,2,4,8,  and  16. 
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Histogram  frequencies  were  found  and  compared  with  the  frequency 
distribution  for  random  normal  deviates.  There  appeared  to  be  no 
significant  discrepancy  between  the  histogram  and  the  theoretical 
curve  for  In  all  work  reported  here  M*  10  was  used. 

To  simulate  errors  having  a  given  covariance  matrix  S,  we  first 
determine  [  17]  a  matrix  W  such  that 

S  =  W  rW  .  (40) 

T 

If  V  =  (Vr  V  ,  ...  V  )  is  a  sample  of  independent  random  normal 
~  T 

deviates,  we  form  V  =  W  V.  Then,  since  S^.  =  I,  it  follows  from  the 
propagation  rule  for  covariance  matrices  that 

S~  =  (WT)  I(WT)T  =  WTW  =  S.  (41) 

This  process  was  used  to  simulate  errors  having  zero  mean  and 
covariance  matrix  S^. 

In  the  simulation  of  errors  the  mean  vector  f  was  varied  rather 

o 

than  the  sample  vector  f  because  the  integration  to  give  the  right 
side  of  (36)  could  then  be  obtained  in  closed  form  as  (38),  and 
a  check  on  the  accuracy  of  quadrature  could  also  be  obtained.  A 
total  of  15  quadrature  abscissas  were  chosen,  these  consisted  of  5 
Gass -Legendre  values  for  each  of  3  intervals:  0  to  1,  ltoS^and  3  to  10 
respectively.  Intervals  of  different  length  were  needed  because  of 
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the  general  decay ing -exponential  character  of  the  kernel.  The 

interval  enclosing  the  x.  was  taken  as  1  Sx,  ^  1 . 98,  i  =  1 .  .  .  n.  The 

i  i 

x^  were  evenly  spaced,  and  the  values  used  will  be  indicated.  Under 

these  conditions  a  numerical  integration  with  Gaugsian  quadrature  was 

compared  (for  f-f  )  with  (38)  and  the  maximum  quadrature  error 

was  found  to  be  .0011  in  the  right  side,  for  which  the  actual  value 

varied  from  about  103  to  292.  All  solutions  used  the  classical 

Obukhov  basis.  In  table  2  results  are  shown  for  x.  =  1(.0£44...)  1.98, 

where  the  x^  are  rounded  to  four  decimal  places.  Here  n=  19,cl(x.) 

2  -4 

is  given  by  (37),  is  as  given  in  table  1,  S  =o„  I- 10  -I 

x  €  t  ^ 

and  f  =293.997  as  indicated.  The  solutions  were  run  withk-5. 
c 

Under  these  conditions  TrSTI_  =  TrX  *  =  12.8.  There  was  some 
loss  of  accuracy  in  the  computation,  so  that  trace  values  are  good 
only  to  one  decimal  place.  These  were  the  only  trace  computations 
for  which  the  computation  errors  were  noticeable;  most  other 
trace  values  were  computed  vising  double  precision.  The  overall 
rms  error  corresponding  to  TrX  *  =  12.8  is  .92.  The  solution  of 
table  2  and  two  others  are  shown  graphically  in  figures  7,  8,  and  9. 
The  three  solutions  were  selected  at  random  from  a  "run"  of  40 


solutions. 


Table  2 


SOLUTION  HI  CF  INTEGRAL  EQUATION,  f  =  293.997; 

c 

SAMPLE  RMS  ERROR:  0.71 


y 

UC 

£  =  £  AUC 
o 

£ 

o 

.  046910077 

-1. 220 

294. 02 

295.24 

.23076534 

O' 

CO 

H 

I 

293.  25 

295. 15 

. SOOCOOOO 

-2. 687 

293. 86 

296. 55 

. 7  69234o6 

-3. 018 

295.  07 

298.09 

. 95308992 

-3.  C75 

295. 03 

298. 10 

1. 0938202 

-3. 066 

294. 89 

297. 95 

1.4615307 

-2. 944 

294.05 

297. 00 

2.0000000 

-2.475 

293, 85 

296. 33 

2.5384693 

-2. 155 

293. 60 

295.75 

2.9061798 

-1. 344 

293.44 

295. 28 

3.  3283705 

-1. 59S 

293.98 

293. 58 

4. 6153574 

-1. 157 

292.94 

294. 10 

6. 5000000 

-1.779 

295. 16 

296. 94 

8. 3846426 

-2. 332 

293. 55 

295. 89 

9. 6716295 

-2.268 

294. 52 

296.79 
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The  index  TrX  was  computed  for  exactly  the  sane  conditions  as  for 

the  solutions  of  figures  7,8,  and  9  except  that  S  -■  ^  -  I=(.01}-I-  Under 

these  conditions  it  was  found  that  TrX  1  =  21.44.  Next  every  second  and 

third  value  of  x.  was  removed  to  give  the  set  x.  =  1  (.  16?3.  .  . )  1. 98 
l  i 

where  the  resulting  x,  are  rounded  to  four  decimals.  In  this  case 

°  l 

n  =  7  and  TrX  *  =  23.93.  This  is  the  case  shown  in  figure  I. 

Decreasing  n  from  19  to  7  had  very  little  effect  on  the  expected  error. 

The  effect  of  changing  n  is  summarized  in  table  3.  For  this  kernel,, 
the  theoretical  rms  error  is  not  strongly  influenced  by  n,  as  long  as 
the  interval  (1,  1.98)  remains  fixed. 


Table  3 

TrX*1  FOR  VARIOUS  x -CONFIGURATIONS  ON  (1,  1.  98) 

n  x.  TrX  1  Theoretical 

Overall  rms  errors 


19 

1(.  0544)1. 98 

21.44 

1.20 

7 

1(. 1633)1. 98 

23.  95 

1.26 

2 

1(.  98)  1.98 

27.45 

1. 35 

Figures  10,  11,  12,  and  13  present  results  for  the  conditions  of  tabic  3 
with  n  =  7  and  k  =  4.  These  solutions  were  selected  at  random  from 


a  "run"  of  40  solutions.  The  improvement  of  the  solution  over  the 
statistical  values  is  not  as  pronounced,  because  of  the  increase  of  <r { 
from  .01  to  .1. 
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Table  4 

TYPICAL  SOLUTION  OF  INTEGRAL  EQUATION  FOR 
VARIOUS  VALUES  OF  k  WITH  n  =  7  and  0  =  0.  1 


The  40  solutions  of  which  figures  10  through  13  are  a  sample 
were  also  run  for  k  =  1 ,  2,  3,  5,  6,  9»  1  2,  and  15.  The  same  sequence 
of  random  errors  was  used  for  each  run  of  40.  A  typical  solution  is 
shown  in  table  4  for  various  values  of  k.  Table  4  shows  typical 
behavior  of  the  solution  as  k  is  varied.  In  this  case,  any  value  of 
k  from  3  through  15  would  have  given  essentially  the  same  solution. 


7.  Absorption  of  Microwaves  in  the  10-150  GHz  Frequency  Band. 


The  Van  Vleck-Weisskopf  theory  [  21]  was  applied  by  Van  Vleck 
[22]  ,  f  2 3]  ,  to  the  absorption  of  microwave  by  oxygen  ana  water  vapor. 
The  formula  obtained  by  Van  Vleck  [22]  for  the  absorption  of  is 

,  8tt3  v  N  ,  Rc  B 

y  =  10  (log ^ 0  c) - — —  2 


c(3k  T) 

Lf2|W‘Wv|lbJ2  +  2|V 


N 


N- 


+  F(v)':  pNl  ]  e 


-E 


N/kT 


(42) 


where  the  following  notation  has  been  introduced: 

a)  is  the  total  number  of  molecules  per  cubic  centimeter 

and  R  is  the  fraction  of  the  total  which  is  0_; 

°2  2 

b)  c,  k,  and  v  are  the  speed  of  light,  Boltzmann's  constant,  and 
the  frequency  of  the  imposed  microwave,  respectively; 

c)  f  and  F  are  form  factors  for  the  line  shapes  given  below; 

d)  VN+  and  v  are  resonant  frequencies  of  the  transitions; 

(N0  labels  the  so  called  zero  frequency  transition)* 

e)  V  is  the  absorption  in  decibels  per  kilometer  and  is  numerically 
equal  to  the  product  of  tne  mass  absorption  coefficient  times  the  den 

f)  N  runs  over  odd  integers  labeling  certain  rotational  quantum  states; 

2 

{?)  the  square  of  tne  matrix  clement  dipole  moment  evaluated 

between  appropriate  angular  momentum  states; 
h)  F>  is  a  relational  constant  whose  magnitude  is  1.44  cm  \  and 


The  u  amnion  frequencies  for  oxygen  can  be  calculated  with  sufficient 
accuracy  for  tne  purposes  of  this  report  from  the  formulas  in  Townes  and 
Schawlow 
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i)  £n  is  the  energy  level  of  the  N*h  quantum  state. 

In  (42)  the  expression  of  Van  Vleck  has  been  translated  into  the 
more  modern  notation  and  two  of  these  equations  are  consolidated. 
Furthermore, VanVleck’s  expression  is  for  pure  oxygen.  The  fraction, 
Ro2,  of  oxygen  in  the  dry  atmosphere  is  about  0.2.  In  the  actual  atmos¬ 
phere,  R  varies  slightly  since  the  mixing  ratio  of  water  vapor  also 
2 

varies. 

The  form  factors  in  (42)  are  given  by 


£(v.,v)  =  __[ 


i  (V.-V)  +  (2sv)  {v  +v)6  +  (4V)‘ 


lim  [vf(v,v)J 
v.-o 

i 


2  v  Av _ 

ir  [  V2  +  (£v)2  ] 


The  matrix  elements  of  the  dipole  moment  can  be  calculated  from  form 
uias  in  Townes  and  Schawlow  [  24], 

The  formula  for  the  absorption  by  water  vapor  at  the  22.235  GHz 
resonant  line  in  decibels  per  kilometer  is  given  by  Van  Vleck  (  23]  as 

v,  „  6  8’3vNd8H20  2  -Vl/kT 

6.5  5_j  =  10  (log10e)  y-g-jjp -  5  I  e 


and  G  is  the  rotational  partition  function  given  by 


P  J  ,  ,  "EJT/kT 

°  =  I  I  [2  ~  (-1)  T  ]  (2J+1)  e 


J=0  T= -J 
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The  dipole  moments  and  transition  frequencies  of  H^O  may  be  calcu¬ 
lated  from  formulas  in  the  Appendices  of  Strandberg  [25].  These  are 
much  more  complicated  than  the  corresponding  calculations  of  Cb,. 

Equation  ^46)  implies  that  other  transitions  arc  possible  for  the 
H^O  molecule.  If  the  form  factor  in  (43)  were  valid  for  these,  then 
the  absorption  of  each  of  these  lines  would  have  the  same  form  as  (45). 

These  frequencies  arc  much  higher,  ;  ,  general,  than  the  highest  frequency 
of  interest  for  this  paper.  Thus,  an  approximate  formula  for  these  con¬ 
tributions  Y  can  be  written  as  ,  given  by 

Y„  =  K  p  v2  \47 ) 

Res. 

The  total  Y  due  to  H  O  is  the  sum  of  the  contributions  of  (45)  and 
£ 

(47). 

The  experimental  values  appropriate  for  the  K  and  Av's  for  the  Van  Vleck 
equations  will  be  discussed  in  the  next  section.  The  Av  is  derived  in  the 
VanVleck  -  Weis  skopf  theory  from  considering  broadening  to  arise  from  the 
finite  extension  of  the  time  between  collisions.  In  fact,  2n£v  is  the  recip¬ 
rocal  of  the  mean  time  between  collisions  in  the  VanVleck-W eisskopf  theory. 
If  one  had  hard  spherical  molecules,  the  mean  time  between  collisions 

could  easily  be  calculated  from  kinetic  theory  l  26]  .  However, 
molecules  interact  in  a  much  more  complicated  way  than  can  be  described 

by  the  simple  hard-sphere  collision  picture.  A  better  approximation  can 
be  made  bv  considering  the  molecular  interactions  to  be  more  compli¬ 
cated  functions  of  the  separation  than  the  hard-sphere  approxi¬ 
mation.  This  at  least  improves  the  temperature  dependence  of 


34 


The  classical  theory  described  above  does  not  explain  why  different 
lines  have  different  widths.  One  could  account  for  this  in  the  classical 
approximation  by  ascribing  a  different  interaction  between  each  state. 

No  a  priori  knowledge  from  classical  theory  can  supply  this  sort  of 
information  and  it  must  be  introduced  as  an  additional  set  of  assumptions 
if  one  is  to  proceed  further  via  a  classical  mechanical  calculation  of 
line  widths. 

A  partial  quantum  mechanical  approach  to  the  problem  was  given 
by  Anderson  as  quoted  in  Tsao  and  Gurnutte,  [  27],  The  translational 
motion  in  this  calculation  is  treated  as  a  classical  mechanical  problem, 
whereas  the  internal  state  transitions  are  treated  as  a  quantum  mechani¬ 
cal  problem.  In  this  picture  A  v  is  given  by 


where 


Av  =  n  v  0  , 


a 


b  db  f(I,  b)  . 


(48) 


In  (48)}  n  is  the  number  density  of  perturbers,  v  is  the  average  relative 
velocity,  b  is  the  separation  of  the  molecules  at  closest  approach  and 
f(I,  b)  is  the  distribution  function  for  cross  section  as  a  function  of  the 
internal  coordinates,  I,  and  b. 
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8.  Least-Squares  Analysis  of  Line  Width  Data 

The  absorption  equations  given  in  the  last  section,  as  developed  by 
VanVleck  and  Weisskopf  [21],  contain  essentially  one  undetermined 
parameter,  the  line  width  constant,  A  v/c .  Calculations  of  this  quantity 
are  laborious  and  many  uncertainties  are  present  in  such  calcul«-..ons. 
Examples  of  such  uncertainties  are  lack  of  knowledge  of  the  exact  form 
of  the  molecular  interaction,  imprecise  knowledge  of  molecular  constants 
(such  as  quadrapole  moments)  which  enter  into  the  theory  of  line  broad¬ 
ening,  and  the  effect  of  higher  order  corrections  to  the  perturbation 
expansions. 

As  a  check  on  theory,  there  are  fortunately  many  experimental 
results  on  line  broadening  and  absorption  coefficients  for  a  variety  of 
conditions  of  pressure  and  temperature.  For  low  pressures  (around 
10  to  20  mm  of  Hg),  direct  measurements  of  line  widths  can  be  made 
since  unresolved  lines  are  well  separated  and  there  is  very  little  line 
shape  asymmetry  from  (he  "negative  frequency"  portion  of  the  VanVleck- 
Weisskopf  line  shape.  However,  there  is  danger  in  extrapolating  these 
low  pressure  measurements  to  high  pressures,  especially  if  the  individ¬ 
ual  lines  are  not  well  resolved  at  the  higher  pressures.  A  notable 
example  of  this  type  of  behavior  is  oxygen. 

In  the  work  reported  in  this  section,  the  absorption  data  of  various 
researchers  have  been  least-square  analyzed  to  determine  line  widths 
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and,  in  the  case  of  water  vapor,  line  shape  factors.  A  least- square s 
analysis  applied  to  nonlinear  equations  (such  as  the  VanVleck- Weisskopf 
equations)  cannot,  in  general,  be  solved  exactly;  the  solution  must  be 
determined  iteratively  by  use  of  a  reasonable  first  estimate  to  insure 
convergence. 

The  least-squares  method  used  here  may  be  described  as  follows. 


We  wish  to  characterize  N  data  points  d.,  by  a  function  f.  of  n  parameters 


p.  such  that 
J 


y  (d.  -  f.(p,.p,»  •  •  *P  ))  w.  •-  min  , 
Lj  i  i  i  d.  n  i 

i=l 


where  w  =  positive  definite  weight  of  tjye  i  data  point.  It  is  assumed 

that  "reasonable"  first  approximations  to  the  parameters  p.  are  known. 

If  these  first  approximations  are  labeled  as  p.^  ,  then  the  functions  f 

J  i 

may  be  expanded  in  a  Taylor  series  about  p.(°>  . 

W  ‘  VP'0’»  ♦  l  t  -  p/0’)  • 


5f! 

3h 

3P1 

3P 2  '  '  ' 

3p„ 

Sf2 

3f2 

3f2 

°P1 

3Pz  '  " 

3pn 

?fN 

•fN 

5  f 

N 

Spl 

"N  *  *  • 

=  P2 

- n - 

5pn 
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and 


r 


Ui 


V  d2 


2’ 


,dN  '  fN-* 


(52) 


In  (51)  and  (52), 
evaluated  at  the 


=  vector  of  residuals. 


both  the  function  f.  and 

l 

original  estimate  p. 


its  derivatives  are  to  be 
Then  minimizing  (49)  yields 


the  matrix  equation  given  by  Strand  [  28]. 


ap  P  -  P(0)  -  uTr  j)'1  jTrr> 


(53) 


where 


T  =  diag  (wr  wz>  ....  wN)  .  (54) 

Equation  (53)  is  then  iterated  by  reevaluating  (51)  and  (52)  at  each 
stage  of  the  iteration  process- 

Equation  (53)  was  applied  to  the  oxygen  absorption  data  of  Artman 
[  29]  and  Crawford  and  Hogg  [  30]  to  determine  the  values  of 
the  line  width  constants  for  both  the  (+)  and  (-)  transitions.  The  data 
of  Artman  were  taken  under  laboratory  conditions  at  T  =  300° K  and 
pressures  of  1  atm,  ?  atm,  and  4  atm  in  air.  Each  of  the  above  samples  of 
data  was  taken  at  .approximately  25  frequencies  covering  the  range  of 
50  to  62  GHz. 
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Each  piece;  of  data  was  weighted  with  unit  weight.  The  data  of  Crawford 
and  Hogg  wore  taken  from^transmission  studies  through  the  atmosphere 
at  sea  level  pressures.  As  stated  in  their  article,  the  approximate 
temperature  range  during  the  period  of  measurement  was  0  ±  3°C. 

Hence,  as  with  most  data  taken  in  the  field,  the  uncertainty  of  these  data  is 
greater  than  tnat  of  laboratory  measurements.  The  results  of  the  least-squares 
fit  are  shown  in  table  5.  As  Tftts  been  discussed  previously  [31], 
the  numerical  values  of  the  normalized  line  width  parameter  (A  v/c)/P 
decrease  with  increasing  pressure.  A  new  feature,  possibly  of  significance 
for  line  breadth  theory,  is  that  the  line  width  §  of  tnc  +  transitions  are 
consistently  higher  than  those  of  the  -  transitions  for  all  of  the  data 
examined.  In  addition,  the  data  of  Crawford  and  Hogg  give  higher  values 
of  the  line  width  than  a  simple  extrapolation  using  the  usually  assumed 
temperature  dependence  would' indicate.  This  fact  causes  the  authors  to 
speculate  that  the  interaction  of  H  O  with  might  give  rise  to  the 
increase  in  oxygen  line  widths.  Of  course,  the  residual  absorption  of 
the  22  GHz  and  the  183  GHz  H^O  lines  might  also  appear  in  this 
analysis  as  an  increase  in  line  width,  but  because  of  the  absence  of 
relative  humidity  measurements  while  the  absorption  measurements  were 
being  taken,  it  was  impossible  to  ascertain  if  this  was  the  case.  In  any 
case,  it  would  be  useful  if  measurements  were  taken  of  oxygen  line 
broadening  in  the  presence  of  water  vapor. 
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The  22 -GHz  water  vapor  absorption  data  of  Becker  and  Autler  [32] 
were  also  analyzed,  but  in  a  slightly  different  fashion.  For  a  fixed  temp¬ 
erature,  the  water  vapor  absorption  coefficient  may  be  written  [23]: 


=  C.V  { 

p  1 


Av 


Av/c 


(i  -  v0)2  +  (Av/c)2  (v  +  vor+  [—) 


_ )  +  c  V  — , 

2  ,Av,2  J  2  c 


(55) 


where  — —  =  a(l  +  b  p).  From  66  measurements  of  vs.  p  and  v  ,  four, 
c  p 

and  five-parameter  fits  to  the  data  were  determined.  The  parameters 
determined  were  c.,  c  ,  a,  b,  and,  in  addition,  for  the  five-parameter 

A  Cm 

>. 

fit  VQ.  The  results  are  tabulated  in  table  6  .  The  initial  estimates  in  all 
cases  were  taken  from  accepted  values  of  these  constants. 

Note  that  the  resonant  frequency  determined  by  the  five-parameter  fit  - 
is  very  close  to  the  accepted  value  .7417  cm  *,  which  is  determined 
by  other  means.  The  values  of  c  and  b,  which  describe  respectively 
the  nonresonant  contribution  of  the  183  GHz  and  higher  lines  and  the 
effect  ofH  O-H  O  line  broadening,  are  both  notic eably  higher  than  the 

Ct  Cm 

accepted  values.  Again,  these  facts  could  be  of  theoretical  significance. 

Surprisingly  good  consistency  was  obtained  between  the  four-and-five- 
parameter  determinations,  considering  that  the  data  were  somewhat 
noisy. 
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Table  6 


Water  Vapor  Absorption  Parameters  Determined  by  a  Least-Squares  Analysi 
of  the  Data  of  Becker  and  Autler  Using  the  Form 

av/c _  Av/c  i  ~2  Av 


?  '  V* 


(v-v/  4  (AV/c)2 


_ Av/c _ 

(v+vo)2  4  (Av/c) 


4c2V 


where  Av/c  =  a(l  4  b  p) 


Units:  y/p  dB/km  per  m ' 

/  -1 
Av/c  cm 

-1 

a  cm 


(g/m3)  1 


-1 


cm 


-Parameter 

Fit 

Initial  Estimate 

10th  Iteration 

a 

.  087 

.08478 

b 

.005 

. 00708 

ci 

.004 

.00361 

C2 

.  05 

.06089 

rms 

average 

deviation  in  y/p  =  . 
y/o  =  .0154  dB/k 

00114  dB/km  per 
m  per  g/rri 

i-Parameter 

Fit 

Initial  Estimate 

,  „  th  _ 

10  Iteration 

a 

.087 

. 08366 

b 

.005 

. 00710 

C1 

.  004 

. 00355 

C2 

.05 

. 06204 

vo 

.  7417 

. 74380 

rms  deviation  in  y/p  =  .00114  dB/km  per  g/m' 
average  y/p  =  .  0154  dB/km  per  g/m^ 
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9.  Application  of  Quality  Criterion  to  Determine  the  Reduction 
of  the  Variance  of  a  priori  Temperature  Data  by 
Microwave  Radiation  Measurements 

In  recent  work  by  Twomey  [  ll]  and  Mateer  [8],  the  information 
content  of  radiation  measurements  has  been  discussed.  In  particular, 
it  was  shown  that  for  physically  realistic  smooth  kernels,  independent 
information  obtained  from  separate  measurements  is  severely  limited 
due  to  ubiquitous  measurement  error.  Twomey  [ll]  goes  further  in 
his  analysis  and  determines,  from  a  large  number,  p,  of  possible 
measurements  with  a  given  noise  level,  the  number  of  independent 
pieces  of  information,  £  ,  which  can  be  found  [  ll],  and  gives  a  pro¬ 
cedure  for  determining  optimum  frequency  locations  for  the  ^measure¬ 
ments  (£  < p).  In  many  inversion  problems,  however,  certain  informa¬ 
tion  about  the  unknown  function  is  known  before  any  measurements  are 

made.  This  a  priori  knowledge  is  embodied  in  the  mean  and  the 
covariance  matrix,  both  of  which  can  be  estimated  from  past  data. 

The  information  content  of  radiation  measurements  should  be  judged  by 
the  improvement  over  the  a  priori  knowledge. 

As  noted  in  section  5,  the  quantity  TrX  1  can  be  used  as  a 

measure  of  the  overall  accuracy  to  be  expected  from  a  given  integral 
equation  with  given  covariance  matrices  Sf  and  S  and  a  given  set 
ix^i  of  x.  Calculating  TrX  ^  and  comparing  with  TrS^  indicate 
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the  amount  of  improvement  over  the- a  priori  statistical  knowledge 
that  can  be  expected  in  a  given  case.  It  is  important  to  realize, 
however,  that  TrX  is  not  a  panacea  because  no  single  number 
can  adequately  describe  m(m  +  l)/2  of  them.  More  detailed  informa- 

tion,  such  as  the  amount  of  reduction  of  variance  of  the  function  at 

^  -1 
various  points,  can  only  be  obtained  by  comparison  of  and  X 


The  calculations  of  X  in  this  section  have  been  carried  out  for 
ground-based  probing  "looking  up"  and  for  an  integral  equation 

V 

determined  by  the  Van  Vleck  equation  (42)  for  microwave  oxygen 
absorption.  The  temperature  and  pressure  structure  used  in  calcu¬ 
lating  the  absorption  at  each  height  were  determined  by  assuming  a 
constant  lapse  rate  of  6.5  °K/kmto  determine  the  temperature  and 
by  integrating  the  hydrostatic  equation  to  determine  the  pressure. 

The  effect  of  non-linearities  in  the  kernel  due  to  fluctuations  of  the 
temperature  from  the  mean  value  has  not  been  included  in 
this  analysis.  The  water  vapor  absorption  was  assumed  to  be  zero, 
and  surface  temperature  and  pressure  were  taken  as  0°  C  and  850  mb, 
respectively.  For  simplicity  the  pressure  dependence  of  the  normalized 
line  width  constant  (Av/P)  was  neglected  and  the  values  of 
^(-£-)/PJ±  =  029g  cm  /atm  were  used  for  +  and  -  transitions  respec¬ 

tively.  The  calculations  at  each  frequency  have  been  made  assuming 
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m  onoch  r  ora  at  ic  ity ,  so  that  the  degradation  oi  information  due  to 
finite  receiver  bandwidth  is  not  accounted  for.  This  effect  should 
be  small,  however,  since  the  ratio  of  receiver  bandwidth  to  pressure- 
broadened  linewidth  can  be  made  quite  small  for  microwave  receivers. 
Because  of  the  previous  assumptions,  the  calculations  involving  X  ^are 
meant  to  be  illustrative  only. 

Before  presenting  these  calculations,  we  must  give  one  more  theo¬ 
retical  extension  of  the  methods  of  sections  3,4,  and  5.  For  a  ground- 
based  probing  scheme,  the  value  of  the  function  at  the  lower  point,  a, 
can  usually  be  measured  directly.  This  constrained  point  should  be 
used  to  modify  the  statistical  estimation  of  the  function  and  its 
uncertainty  as  given  by  X  \  First,  the  constraint  can  be  incorporated 
into  the  integral  equation  by  using  a  quadrature  formula  (such  as 
Gauss -F.adau)  which  uses  the  value  of  the  function  at  a,  f(a),  directly: 

b  m 

C  K(x,  y)  f  (y)  dy  te  V  w  K(x, y  )  f(y.) 

j  =  1  J  J  J 

m 

=  Wj  K(x ,  a)  f  (a)  +  Y  w.  K(x,  y.)f(y.) 

i  -  ?  J  J  ) 

j  *• 

Thus,  by  subtracting  w^  K(x  ,  a)  f  (a)  from  the  measured  quantity, 
g(x),  a  matrix  equation  to  be  solved  for  (m-1)  components  of  the 


(56) 


45 


'  ■  M"*! OMLa. 


function  of  f  is  obtained.  Second,  knowing  f{a)  redut  es  the  uncer- 

tainty  in  all  the  other  functional  values.  The  new  covariance  matrix, 

(c ) 

S ^  ,  of  dimension  (  m  -  1)  x  (m  -  1),  has  elements 

lj  j  “  J  ^ |  •  •  •  j  9 


s(c)  =  S 

ij  M 


S,  .  S.. 

JLlJLl 


'n 


where 


S.  =  (S.  )  and  S,(c)  -  (S..(C))  . 
i  y  f  ij 


(57) 


For  convenience,  the  matiix  S 


(c> 


f 


will  be  referred  to  in  the  follow¬ 


ing  as  the  constrained  covariance  matrix.  Furthermore,  instead  of  the 
mOn,  f,  as  the  best  a_  priori  estimate  of  f,  the  effect  of  knowing  f^ 
modifies  this  optimum  a  priori  estimate  to  ?  where 


-  -  Sil 

f.  =  f  .  +  -±-  (f  -  f  ),  i=  2,  ....  m. 

i  1  1  1 


(58) 


Equations  (57)  and  (58)  may  be  derived  from  linear  regressions  of  the 

(m  -  1)  functional  values  f  ,  f f  as  functions  of  the  surface  value 

2  3  m 

£r 

The  two  temperature  covariance  matrices  and  the  corresponding 
means  were  obtained  by  ensemble  averaging  of  5  years  of  radiosonde 
data  for  the  ?,  months,  February  and  August,  at  Denver,  Colorado, 

The  February  and  August  averages  were  over  163  and  240  soundings, 
respectively.  The  temperatures,  given  at  heights  determined  by  the 
radiosondes,  were  interpolated  linearly  to  determine  values  at  the  15 
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desired  quadrature  points.  (Gauss -Radau  quadrature  was  used  for  the 
subinterval  0  to  1  km,  while  Gaussian  quadrature  was  used  in  the  sub¬ 
intervals  1  to  3  km  and  3  to  10  km  .)  The  means  T  (h^),  and  elements  of  the 
covariance  matrix,  S(h.,  h.)  s  s...  were  determined  fr< 

i  j'  ij’ 


•om 


and 


^  N 

*  “V =  i  I  <>h)  ■ 


(59) 


% 


N 

ibi  1 i  (T° 


(60) 


where  N  is  the  number  of  soundings,  a  is  an  index  for  each  sound¬ 
ing,  and  Ta  (hj)  is  the  temperature  of  the  sounding  at  the  height 
h^.  The  constrained  covariance  matrix  was  then  calculated  using 
(57).  The  means  of  February  and  August  are  .given  in  table  7  and 
the  constrained  covariance  matrices  for  February  and  August  are 
given  in  tables  8  and  9,  respectively.  It  may  be  of  interest  to 
compare  the  unconstrained  covariance  matrix  for  August,  given  in 
table  1,  with  the  constrained  matrix  for  the  same  month,  table  9. 
The  elements  of  the  matrices  differ  considerably  only  below  the 
10^  quadrature  point  (around  3  km). 
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Tabic  7 


Five-year  Mean  Temperature  Profiles  at  Denver,  Colorado, for  February 
and  August . 


i 

hr 

_  (km) 

Tlh.)  (Feb.  ) 

(°K) 

T  (h^  (Aug.) 
(°K) 

i 

0.  000 

268. 23 

' 

294. 00 

2 

0.  140 

269.  92 

293.48 

3 

0.416 

269.  56 

292. 47 

4 

0.723 

268.  15 

290. 87 

5 

0.  943 

266.  83 

289.43 

6 

1.094 

265.  87 

288.40 

7 

1.462 

263.  54 

285. 69 

8 

2.000 

260. 08 

281. 56 

9 

2.  538 

256.76 

277.32 

10 

2.906 

254.40 

274.47 

11 

3.  328 

251. 59 

271.25 

12 

4.  615 

242. 56 

262.74 

13 

6.  500 

228. 80 

250. 11 

14 

8.  385 

219. 18 

236. 16 

15 

9.  672 

217.11 

226.76 

To  determine  the  amount  of  reduction  of  the  a  priori  temperature 
variance  by  microwave  measurement;  of  oxygen  emission,  calcula¬ 
tions  of  X  1  were  made  for  a  set  of  frequencies  which  approximately 
cover  the  entire  oxygen  band.  This  set  is  shov/n  in  table  10. 


Temperature  Covariance  Matrix,  in  (°K)  ,  of  Denver,  February,  Radiosonde  Data  for 
the  2nd  to  15t’n  Quadrature  Points  with  Constrained  Surface  Temperature. 
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Table  9 

Temperature  Covariance  Matrix,  in  (  K)  ,  of  Denver,  August,  Radiosonde  Data  for 
the  2nd  to  the  15th  Quadrature  Points  with  Constrained  Surface  Temperature. 
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Table  10 

Frequencies  Used  in  Calculations  of  T  r  X  ^ 
v(GHz) 

47.02650 
47.22650 
47. 94917 
48.45304 
50.28294 


52. 02593 
53. 93117 
55. 22163 
56.  26466 
58.44669 


60.43505 
61. 80036 
62. 48631 
62.  68631 
63. 98631 


The  calculations  of  X  were  run  for  both  months,  and  for  assumed 
rms  brightness  temperature  measurement  errors  of  .  01  °K,  .  1  °K, 
and  1.0  °K.  The  entire  matrix  of  X  *  for  Denver,  February,  with 
o,=  .01  °K  is  shown  in  table  1 1 . 

The  meaning  of  a  15x15  covariance  matrix  cannot  be  presented  in 
simple  form.  A  rough  estimate  of  the  standard  deviation  to  be 
expected  at  each  quadrature  height  is  given  by  the  square  root 
of  the  corresponding  diagonal  element  of  the  covariance  matrix. 

These  quantities  are  plotted  as  functions  of  height  in  figures  14  through 
19.  For  ease  of  comparison,  the  error  distributions  to  be  expected 
in  the  constrained  solution  are  shown  in  figures  20  and  21.  The  cor¬ 
responding  overall  quality  criterion,  TrX  ,  is  tabulated  in  tables 
12  and  13. 
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Table  1Z 


Table  13 
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The  study  of  X  ^  and  its  trace  can  be  directed  towards  an 
important  experimental  consideration  i.  e.  ;  namely,  that  of  determining 
an  optimum  distribution  of  measurement  trequencies,  given  that  only 
a  certain  number  of  them  (usually  small)  will  be  available.  This  problem 
has  been  studied  by  Twomey  [11  J  from  the  point  of  view  that  the 
optimum  location  of  measurements  is  determined  by  the  integral 
equation  itself.  Furticr,  from  the  analysis  of  a  certain  matrix,  which 
depends  on  the  kernel,  he  shows  that  the  number  of  independent  pieces  of 
information  is  limited  by  the  number  of  eigenvalues  which  are  larger  than 
the  noise  level  of  the  measurements.  In  the  absence  of  additional 
information,  the  above  approacn  seems  reasonable. 

When  additional  information  is  available,  such  as  the  statistical 
information  contained  in  S^.,  it  is  cloar  that,  the  optimum  frequencies  should 
be  determined  from  a  statistical  compromise  among  S^.,  S^,  and  the  proper¬ 
ties  of  A.  This  optimum  set  could  be  different  for  periods  of  different 
meteorological  statistics.  If  the  measuring  instruments )  is  (are)  not 
tuneable,  a  heterogeneous  ensemble  averaging  over  yearly  data  might  be  in 
order.  A  method  which  properly  weights  all  available  information  has  the  ad¬ 
vantage  of  allowing  the  choice  of  frequencies  such  that  the  duplica.fion  of  easily 
available  information  (that  of  and  f  )  by  difficult  experimental  measure¬ 
ments  is  made  as  small  as  possible.  The  concept  of  "number  of 
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independent  pieces  ol  informal iun"  is  not  as  straightforward  with 


this  analysis  as  it  was  in  Twumey's.  However,  an  '  effective  number", 
such  as  the  number  of  measurements  required  to  reduce  the 
a_  priori  variance  to  a  desired  level,  given  that  the  noise  level  is 
fixed,  can  be  defined.  In  a  theoretical  sense,  such  a  number  could 
always  be  obtained  for  any  given  level  greater  than  zero  since  dupli¬ 
cating  measurements  is  equivalent  to  lowering  the  noise  level. 
Practically  speaking,  if  a  reasonably  large  reduction  in  variance 
with  a  few  well-placed  measurements  is  not  possible,  the  experiment 
is  probably  impractical. 

An  attempt  was  made  to  determine  the  approximate  sensitivity 
of  TrX  ^  to  choice  of  and  number  of  frequencies  for  one,  two^and 
three  frequencies  and  for  a  =  .  01°K.  The  one -frequency  search  was 
carried  out  starting  with  47.0  GHz  and  adding  a  .5  GHz  increment 
until  56  GHz  was  reached.  These  resultr,  for  both  the  constrained 
and  unconstrained  data,  are  shown  in  table.  14  for  the  Denver, 

August,  statistics. 

The  frequencies  corresponding  to  the  minimum  trace  are  54.  5  GHz 

-1  7 

for  the  unconstrained  (U)  case  (TrX  =  50.6  (  K)*")  and  4.75  GHz  for 
the  constrained  (C)  case  (TrX  -  37.7  (°K)  ).  The  traces  of 
were  120.  5  (G)  and  68.  5  (U),  (°K)^,  and  the  traces  ol  X  ^  for  the 
15  frequencies  of  table  10  were  5.7(  K)^(U)  and  4.  8(  K)  (C).  Thus 
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Tabic  14 


Trace  X  vs.  Frequency,  Denver,  August  r>  ^  =  .01  °K 
(Single  Frequency,  n~  1)  Trace  =  120.48 


V 

(GHi) 

Tr  X  " 1  (  K) 
Unconstrained 

Tr  X  (  K) 
Constrained 

47. 5 

60.  3 

37.7 

48.  0 

56.  9 

39.  2 

48.  5 

55.4 

39.9 

49.  0 

54. 8 

40.2 

49.  5 

54.  5 

40.4 

50.  0 

54.  2 

40.  6 

50.  5 

53.  9 

40.7 

51. 0 

*3.  6 

40.  9 

51.5 

53.  2 

41.2 

52.  0 

52.  8 

41. 5 

52.  5 

52.  3 

41.9 

53. 0 

51.8 

42.4 

53.  5 

51.2 

43.  1 

54.0 

50.7 

44.  2 

54.  5 

50.  6 

45.  6 

55.  0 

51.  0 

47.0 

55.  5 

51.8 

48.  3 

5b.  0 

52.7 

49.4 

5b 


the  use  of  a  sink’ll,  frequency,  with  this  choice  oi  0  ,  reduces  the 

variance  by  a  factor  of  Z  (roughly). 


The  two-frequency  calculations  were  done  by  fixing  V  -  qy  qjj,, 
and  varying  the  second  frequency,  v  ,  from  48.  3  Gllx  to  60.  0  GHz  in 


1.3  GHz  increments.  These  results  are  shown  in  table  15,  which 


indicates  that  the  variance  obtained  from  the  one -frequency  calculations 
is  roughly  halved  by  the  addition  of  the  second  frequency. 

The  three -frequency  calculations  were  done  by  fixing  =  47  GHz 
and  -  50.9  GHz  and  varying  the  third  frequency,  v^,  from  47.5  to 
58.0  in  .5  GHz  increments.  These  results,  again  for  Denver,  August, 
and  o  =  .01  °K,  are  shown  in  table  16.  It  should  be  observed  that 
TrX  *  is  approaching  that  of  the  15 -frequency  calculations,  and 
that  addition  of  more  frequencies  wiLl  yield  diminishing  returns.  The 
relative  minimum  of  TrX  ^  at  47.5  GHz  is  presumably  due  to  the 
nearby  resonant  line  at  47  4465  GHz,  A  rigorous  search  for 
optimum  frequencies  would  necessarily  examine  the  effect  of  each 
resonant  line  on  TrX  . 


57 


Tal.il  e  15 


Two-Frequency  Calculations  of  Tr  X  for  Denver, 
August,  0  -  .01  °K,  v  =  47.0  GHz 


v2 

(GHz) 

T  r  X  ~ 1 

Unc  unstrained 

Tr  X*1 

Constrained 

48.  3 

22.  3 

18.  1 

49.6 

21.9 

17.7 

50.  9 

21.6 

17.4 

52.2 

21. 5 

17.4 

53.  5 

21.8 

17.9 

54.  8 

22.  6 

19.  0 

56.  1 

24.  3 

21. 0 

57.4 

26.2 

2  3.3 

58.7 

27.4 

24.  8 

60.  0 

27.  9 

25.  3 
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10.  Summary 


A  method  for  the  numerical  solution  of  a  Fredholm  integral 
equation  of  the  first  kind  has  been  derived  and  illustrated.  The 
method  requires  knowledge  of  the  covariance  matrices  of  the  con¬ 
straint  vector  and  the  measurement  error  vector.  Such  knowledge 
is  frequenctly  available  in  physical  problems  when  it  is  desired  to 
derive  values  of  physical  functions  from  integrals  involving  them. 

If  both  covariance  matrices  and  Se  are  scalar,  the  equations  re¬ 
duce  to  those  of  Twomey  [  6],  where  the  optimum  smoothing  parameter 
V  is  given  by  the  ratio  of  variances  between  the  diagonal  elements 
of  and  S^,  respectively.  The  present  method  automatically 
incorporates  the  optimum  amount  of  smoothing  in  the  sense  of  maxi¬ 
mum-likelihood  estimation. 


The  trace  of  the  error  covariance  ^matrix,  TrS^^.  is  used 

to  estimate  the  precision  of  the  solution.  When  k=  m  (i.e.,when  the 


basis  forms  a  nonsingular  mxn),  STT_  reduces  to 

UC  -  T| 

-1  -1  T  -1  -1  -1 

X  =  (S^  +  A  S£  A)  The  positive  number  TrX  is  related 


to  the  error  to  be  expected  in  the  solution  and  is  used  as  a  quality 
criterion.  A  comparison  of  TrX  ^  with  TrS^  indicates  the  amount  of 
information  contained  in  the  integral  equation  with  observation  errors 
determined  by  .  It  is  evident  that  TrX  *  can  be  used  to  study 
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optimization  of  the  spacing  of  observations.  Such  a  study  could  be 
valuable  for  planning  measuring  systems. 

The  method  of  estimation  and  the  use  of  the  quality  criterion 
TrX  *  were  illustrated  by  studying  the  effect  of  various  sequences 
of  random  errors  on  the  solution  of  a  specific  integral  equation. 

The  effect  of  varying  the  number  and  type  of  basis  vectors  was 
studied  for  this  example.  The  results  apparently  agree  with  theore¬ 
tical  predictions,  although  confidence  tests  were  not  run. 

Calculations  of  X  ^  and  its  trace  were  then  performed  using 
certain  meteorological  statistics  of  DenveY,  Colorado,  and  absorption 
coefficients  determined  from  the  Van  Vleck  absorption  equations  for 
the  microwave  oxygen  complex  with  a  linear  temperature  profile. 

These  calculations  were  performed  for  a  set  of  15  frequencies 
selected  to  cover  the  oxygen  band,  and  for  three  sets  of  assumed 
brightness -temperature  rms  errors  of  .01  °K,  .1  °K,  and  1.0  °K. 

The  square  roots  of  the  diagonal  elements  of  X  were  plotted  as 
functions  of  height  in  an  attempt  to  show  the  accuracy  with  which  various 
parts  of  the  tropospheric  temperature  profile  could  be 
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i 

4. 


inferred  from  ground-based  measurements  of  oxygen  emission. 
Questions  of  optimum  choice  of  frequencies  for  this  probing  were 
discussed  and  illustrative  calculations  relevant  to  this  choice  were 
presented . 
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Stondard  Deviation  in  °K 


Square  Roots  of  Diagonal  Elements  of  Covariance  Matrices  of  Temperature 
Statistics,  S and  Solution,  X  vs.  Height  for  Denver,  August, and  for 
Measurement  Errors,  o  ,  of  .1  °K. 


Curve  Labels 

U:  Unconstrained  Lower  Point 
C:  Constroined  Lower  Point 
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Square  Roots  of  Diagonal  Elements  of  Covariance  Matrices  of  Temperature 
Statistics,  S,  and  Solution,  X-*,  vs.  Height  for  Denver,  February*  and  for 

Measurement  Errors,  o  ,  of  .01  K. 

e 


Fig.  17 


83 


Standard  Deviation  in  °K 


Square  Roots  of  Diagonal  Elements  of  Covariance  Matrices  of  Temperature 
Statistics,  S,,  and  Solution,  X'^f  vs.  Height  for  Denver,  February  ,  and  for 
Measurement  Errors,  o  ,  of  .  1  K. 
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Standard  Deviation  in  °K 


Square  Roots  of  Diagonal  Elements  of  Covariance  Matrices  of  Temperature 
Statistics,  S^,  and  Solution,  X“*,  vs.  Height  for  Denver,  February,  and  for 
Measurement  Errors,  ~  ,  of  1.0  '  K. 
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